Expanding the SnoRNA Interaction Network: Conservation of Guiding Function in Vertebrates

Small nucleolar RNAs (snoRNAs) are one of the most abundant and evolutionary ancient group of small non-coding RNAs. Their main function is to target chemical modifications of ribosomal RNAs (rRNAs) and small nuclear (snRNAs). They fall into two classes, box C/D snoRNAs and box H/ACA snoRNAs, which are clearly distinguished by conserved sequence motifs and the type of modification that they govern. The box H/ACA snoRNAs are responsible for targeting pseudouridylation sites and the box C/D snoRNAs for directing 2’-O-methylation of ribonucleotides. A subclass that localize to the Cajal bodies, termed scaRNAs, are responsible for methylation and pseudouridylation of snRNAs. In addition an amazing diversity of non-canonical functions of individual snoRNAs arose. The modification patterns in rRNAs and snRNAs are retained during evolution making it even possible to project them from yeast onto human. The stringent conservation of modification sites and the slow evolution of rRNAs and snRNAs contradicts the rapid evolution of snoRNA sequences. Recent studies that incorporate high-throughput sequencing experiments still identify undetected snoRNAs even in well studied organisms as human. The snoRNAbase, which has been the standard database for human snoRNAs has not been updated ince 2006 and misses these new data. Along with the lack of a centralized data collection across species, which incorporates also snoRNA class specific characteristics the need to integrate distributed data from literature and databases into a comprehensive snoRNA set arose. Although several snoRNA studies included pro forma target predictions in individual species and more and more studies focus on non-canonical functions of subclasses a systematic survey on the guiding function and especially functional homologies of snoRNAs was not available. To establish a sound set of snoRNAs a computational snoRNA annotation pipeline, named snoStrip that identifies homologous snoRNAs in related species was employed. For large scale investigation of the snoRNA function, state-of-the-art target pedictions were performed with our software RNAsnoop and PLEXY. Further, a new measure the Interaction Conservation Index (ICI) was developed to evaluate the conservation of snoRNA function. The snoStrip pipeline was applied to vertebrate species, where the genome sequence has been available. In addition, it was used in several ncRNA annotation studies (48 avian, spotted gar) of newly assembled genomes to contribute the snoRNA genes. Detailed target analysis of the new vertebrate snoRNA set revealed that in general functions of homologous snoRNAs are evolutionarily stable, thus, members of the same snoRNA family guide equivalent modifications. The conservation of snoRNA sequences is high at target binding regions while the remaining sequence varies significantly. In addition to elucidating principles of correlated evolution it was possible, with the help of the ICI measure, to assign functions to previously orphan snoRNAs and to associate snoRNAs as partners to known but so far unexplained chemical modifications. As further pattern redundant guiding became apparent. For many modification sites more than one snoRNA encodes the appropriate antisense element (ASE), which could ensure constant modification through snoRNAs that have different expression patterns. Furthermore, predictions of snoRNA functions in conjunction with sequence conservation could identify distant homologies. Due to the high overall entropy of snoRNA sequences, such relationships are hard to detect by means of sequence homology search methods alone. The snoRNA interaction network was further expanded through novel snoRNAs that were detected in data from high-throughput experiments in human and mouse. Through subsequent target analysis the new snoRNAs could immediately explain known modifications that had no appropriate snoRNA guide assigned before. In a further study a full catalog of expressed snoRNAs in human was provided. Beside canonical snoRNAs also recent findings like AluACAs, sno-lncRNAs and extraordinary short SNORD-like transcripts were taken into account. Again the target analysis workflow identified undetected connections between snoRNA guides and modifications. Especially some species/clade specific interactions of SNORD-like genes emerged that seem to act as bona fide snoRNA guides for rRNA and snRNA modifications. For all high confident new snoRNA genes identified during this work official gene names were requested from the HUGO Gene Nomenclature Committee (HGNC) avoiding further naming confusion

Y3IP1, A nucleus-encoded Thylakoid Protein, Cooperates with the Plastid-Encoded YCf3 Protein in Photosystem I Assembly of Tobacco and Arabidopsis

The intricate assembly of photosystem I (PSI), a large multiprotein complex in the thylakoid membrane, depends on auxiliary protein factors. One of the essential assembly factors for PSI is encoded by ycf3 (hypothetical chloroplast reading frame number 3) in the chloroplast genome of algae and higher plants. To identify novel factors involved in PSI assembly, we constructed an epitope-tagged version of ycf3 from tobacco (Nicotiana tabacum) and introduced it into the tobacco chloroplast genome by genetic transformation. Immunoaffinity purification of Ycf3 complexes from the transplastomic plants identified a novel nucleus-encoded thylakoid protein, Y3IP1 (for Ycf3-interacting protein 1), that specifically interacts with the Ycf3 protein. Subsequent reverse genetics analysis of Y3IP1 function in tobacco and Arabidopsis thaliana revealed that knockdown of Y3IP1 leads to a specific deficiency in PSI but does not result in loss of Ycf3. Our data indicate that Y3IP1 represents a novel factor for PSI biogenesis that cooperates with the plastid genome-encoded Ycf3 in the assembly of stable PSI units in the thylakoid membrane

Developmentally Driven Changes in Adipogenesis in Different Fat Depots Are Related to Obesity

Subcutaneous (sc) and visceral (vis) adipose tissue (AT) contribute to the variability in pathophysiological consequences of obesity and adverse fat distribution. To gain insights into the molecular mechanisms distinguishing vis and sc fat, we compared the transcriptome during differentiation of immortalized adipocytes from murine epididymal (epi) and inguinal (ing) AT. RNA was extracted on different days of adipogenesis (−2, 0, 2, 4, 6, 8) and analyzed using ClariomTM D mouse assays (Affymetrix) covering >214,900 transcripts in >66,100 genes. Transcript Time Course Analysis revealed 137 differentially expressed genes. The top genes with most divergent expression dynamics included developmental genes like Alx1, Lhx8, Irx1/2, Hoxc10, Hoxa5/10, and Tbx5/15. According to pathway analysis the majority of the genes were enriched in pathways related to AT development. Finally, in paired samples of human vis and sc AT (N = 63), several of these genes exhibited depot-specific variability in expression which correlated closely with body mass index and/or waist-to-hip ratio. In conclusion, intrinsically programmed differences in gene expression patterns during adipogenesis suggest that fat depot specific regulation of adipogenesis contributes to individual risk of obesity

Homology-based annotation of non-coding RNAs in the genomes of Schistosoma mansoni and Schistosoma japonicum

<p>Abstract</p> <p>Background</p> <p>Schistosomes are trematode parasites of the phylum Platyhelminthes. They are considered the most important of the human helminth parasites in terms of morbidity and mortality. Draft genome sequences are now available for <it>Schistosoma mansoni </it>and <it>Schistosoma japonicum</it>. Non-coding RNA (ncRNA) plays a crucial role in gene expression regulation, cellular function and defense, homeostasis, and pathogenesis. The genome-wide annotation of ncRNAs is a non-trivial task unless well-annotated genomes of closely related species are already available.</p> <p>Results</p> <p>A homology search for structured ncRNA in the genome of <it>S. mansoni </it>resulted in 23 types of ncRNAs with conserved primary and secondary structure. Among these, we identified rRNA, snRNA, SL RNA, SRP, tRNAs and RNase P, and also possibly MRP and 7SK RNAs. In addition, we confirmed five miRNAs that have recently been reported in <it>S. japonicum </it>and found two additional homologs of known miRNAs. The tRNA complement of <it>S. mansoni </it>is comparable to that of the free-living planarian <it>Schmidtea mediterranea</it>, although for some amino acids differences of more than a factor of two are observed: Leu, Ser, and His are overrepresented, while Cys, Meth, and Ile are underrepresented in <it>S. mansoni</it>. On the other hand, the number of tRNAs in the genome of <it>S. japonicum </it>is reduced by more than a factor of four. Both schistosomes have a complete set of minor spliceosomal snRNAs. Several ncRNAs that are expected to exist in the <it>S. mansoni </it>genome were not found, among them the telomerase RNA, vault RNAs, and Y RNAs.</p> <p>Conclusion</p> <p>The ncRNA sequences and structures presented here represent the most complete dataset of ncRNA from any lophotrochozoan reported so far. This data set provides an important reference for further analysis of the genomes of schistosomes and indeed eukaryotic genomes at large.</p

Phylogenetic distribution of plant snoRNA families

Background: Small nucleolar RNAs (snoRNAs) are one of the most ancient families amongst non-protein-coding RNAs. They are ubiquitous in Archaea and Eukarya but absent in bacteria. Their main function is to target chemical modifications of ribosomal RNAs. They fall into two classes, box C/D snoRNAs and box H/ACA snoRNAs, which are clearly distinguished by conserved sequence motifs and the type of chemical modification that they govern. Similarly to microRNAs, snoRNAs appear in distinct families of homologs that affect homologous targets. In animals, snoRNAs and their evolution have been studied in much detail. In plants, however, their evolution has attracted comparably little attention. Results: In order to chart the phylogenetic distribution of individual snoRNA families in plants, we applied a sophisticated approach for identifying homologs of known plant snoRNAs across the plant kingdom. In response to the relatively fast evolution of snoRNAs, information on conserved sequence boxes, target sequences, and secondary structure is combined to identify additional snoRNAs. We identified 296 families of snoRNAs in 24 species and traced their evolution throughout the plant kingdom. Many of the plant snoRNA families comprise paralogs. We also found that targets are well-conserved for most snoRNA families. Conclusions: The sequence conservation of snoRNAs is sufficient to establish homologies between phyla. The degree of this conservation tapers off, however, between land plants and algae. Plant snoRNAs are frequently organized in highly conserved spatial clusters. As a resource for further investigations we provide carefully curated and annotated alignments for each snoRNA family under investigation

Structured RNAs and synteny regions in the pig genome

BACKGROUND: Annotating mammalian genomes for noncoding RNAs (ncRNAs) is nontrivial since far from all ncRNAs are known and the computational models are resource demanding. Currently, the human genome holds the best mammalian ncRNA annotation, a result of numerous efforts by several groups. However, a more direct strategy is desired for the increasing number of sequenced mammalian genomes of which some, such as the pig, are relevant as disease models and production animals. RESULTS: We present a comprehensive annotation of structured RNAs in the pig genome. Combining sequence and structure similarity search as well as class specific methods, we obtained a conservative set with a total of 3,391 structured RNA loci of which 1,011 and 2,314, respectively, hold strong sequence and structure similarity to structured RNAs in existing databases. The RNA loci cover 139 cis-regulatory element loci, 58 lncRNA loci, 11 conflicts of annotation, and 3,183 ncRNA genes. The ncRNA genes comprise 359 miRNAs, 8 ribozymes, 185 rRNAs, 638 snoRNAs, 1,030 snRNAs, 810 tRNAs and 153 ncRNA genes not belonging to the here fore mentioned classes. When running the pipeline on a local shuffled version of the genome, we obtained no matches at the highest confidence level. Additional analysis of RNA-seq data from a pooled library from 10 different pig tissues added another 165 miRNA loci, yielding an overall annotation of 3,556 structured RNA loci. This annotation represents our best effort at making an automated annotation. To further enhance the reliability, 571 of the 3,556 structured RNAs were manually curated by methods depending on the RNA class while 1,581 were declared as pseudogenes. We further created a multiple alignment of pig against 20 representative vertebrates, from which RNAz predicted 83,859 de novo RNA loci with conserved RNA structures. 528 of the RNAz predictions overlapped with the homology based annotation or novel miRNAs. We further present a substantial synteny analysis which includes 1,004 lineage specific de novo RNA loci and 4 ncRNA loci in the known annotation specific for Laurasiatheria (pig, cow, dolphin, horse, cat, dog, hedgehog). CONCLUSIONS: We have obtained one of the most comprehensive annotations for structured ncRNAs of a mammalian genome, which is likely to play central roles in both health modelling and production. The core annotation is available in Ensembl 70 and the complete annotation is available at http://rth.dk/resources/rnannotator/susscr102/version1.02. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (doi:10.1186/1471-2164-15-459) contains supplementary material, which is available to authorized users

Divergent evolution in the genomes of closely related lacertids, <i>Lacerta viridis</i> and <i>L. bilineata</i>, and implications for speciation

Lacerta viridis and Lacerta bilineata are sister species of European green lizards (eastern and western clades, respectively) that, until recently, were grouped together as the L. viridis complex. Genetic incompatibilities were observed between lacertid populations through crossing experiments, which led to the delineation of two separate species within the L. viridis complex. The population history of these sister species and processes driving divergence are unknown. We constructed the first high-quality de novo genome assemblies for both L. viridis and L. bilineata through Illumina and PacBio sequencing, with annotation support provided from transcriptome sequencing of several tissues. To estimate gene flow between the two species and identify factors involved in reproductive isolation, we studied their evolutionary history, identified genomic rearrangements, detected signatures of selection on non-coding RNA, and on protein-coding genes.Here we show that gene flow was primarily unidirectional from L. bilineata to L. viridis after their split at least 1.15 million years ago. We detected positive selection of the non-coding repertoire; mutations in transcription factors; accumulation of divergence through inversions; selection on genes involved in neural development, reproduction, and behavior, as well as in ultraviolet-response, possibly driven by sexual selection, whose contribution to reproductive isolation between these lacertid species needs to be further evaluated.The combination of short and long sequence reads resulted in one of the most complete lizard genome assemblies. The characterization of a diverse array of genomic features provided valuable insights into the demographic history of divergence among European green lizards, as well as key species differences, some of which are candidates that could have played a role in speciation. In addition, our study generated valuable genomic resources that can be used to address conservation-related issues in lacertids

Overexpressing high levels of human vaspin limits high fat diet-induced obesity and enhances energy expenditure in a transgenic mouse

Adipose tissue inflammation and insulin resistance are hallmarks in the development of metabolic diseases resulting from overweight and obesity, such as type 2 diabetes and non-alcoholic fatty liver disease. In obesity, adipocytes predominantly secrete proinflammatory adipokines that further promote adipose tissue dysfunction with negative effects on local and systemic insulin sensitivity. Expression of the serpin vaspin (SERPINA12) is also increased in obesity and type 2 diabetes, but exhibits compensatory roles in inflammation and insulin resistance. This has in part been demonstrated using vaspin-transgenic mice. We here report a new mouse line (h-vaspinTG) with transgenic expression of human vaspin in adipose tissue that reaches vaspin concentrations three orders of magnitude higher than wild type controls (&gt;200 ng/ml). Phenotyping under chow and high-fat diet conditions included glucose-tolerance tests, measurements of energy expenditure and circulating parameters, adipose tissue and liver histology. Also, ex vivo glucose uptake in isolated adipocytes and skeletal muscle was analyzed in h-vaspinTG and littermate controls. The results confirmed previous findings, revealing a strong reduction in diet-induced weight gain, fat mass, hyperinsulinemia, -glycemia and -cholesterolemia as well as fatty liver. Insulin sensitivity in adipose tissue and muscle was not altered. The h-vaspinTG mice showed increased energy expenditure under high fat diet conditions, that may explain reduced weight gain and overall metabolic improvements. In conclusion, this novel human vaspin-transgenic mouse line will be a valuable research tool to delineate whole-body, tissue- and cell-specific effects of vaspin in health and disease